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Abstract. In this article, we improve the order of the convergence of some finite volume solu- 
tions approximating some second order elliptic problems. 

In one dimensional space, we prove that finite volume approximations of order 0{h''~^^), with k 
integer, can be obtained after k correction using the same scheme of three points and changing 
only the second members of the original system. 

This is done for general smooth second order elliptic problems. These results can be extended 
for non linear second order equation u" = f(x,u,u') where / is a smooth function . 
In two dimensional space, we prove that finite volume approximation of order 0{h?) can be ob- 
tained, starting with finite volume solution of order 0{h), by using the same matrix and changing 
only the second member of the original system. 

This is done for second order elliptic problems of the form —Au + pu = /, with Dirichlet condi- 
tion. 

These results can be extended to obtain finite volume approximation of order 0(h'''^^). 

Heart idea behind these results is the one of Fox's difference correction in the context of finite 

difference method. 
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1. Introduction 

Numerical methods for partial differential equations can be divided into three general categories: 
finite difference methods, finite element methods, finite volume methods. 

Finite difference and finite element methods have been attracted much more attention than finite 
volume methods, consequently there is a well developed literature in finite difference/finite element 
methods which treats several methods for improving the order of the convergence of the approxi- 
mate solutions those using lower scheme. 

The desire to use low order scheme to produce highly accurate approximation in finite difference 
methods led Fox jHj to introduce his difference correction technique. His idea has been modified 
by Pereyra and Lindberg's deferred correction. Theirs ideas have been developed by many authors 
like Zadunaisky's global and Frank's local defect correction (for more informations see j3| and 14 ). 
Almost, the theoretical justifications of these methods are based on the existence of a smooth as- 
ymptotic error expansion for the base scheme. The uniformity of the mesh and the contraction 
property have been the main tools to prove such existence of the error expansions. 



1 



2 



B. ATFEH AND A. BRADJI 



In finite Element methods, defect correction technique has been used to produce highly order of 
convergence by using linear /bilinear finite element method. This has been introduced by Barrett 
et al. and Moore [El in one dimensional space by using uniform mesh and recently by using 
the so-called supraconvergent mesh condition in 

In two dimensional space, under the uniform mesh Chibi [H] ( see also the idea of contraction 
property in this context in Gao et al. jlU| ) has proved that, we can do only one correction on the 
rectangle and corrections we wish for periodic problems ( for a theoretical framework, you can also 
see the communication of Hackbusch in [1], pages 89-113). 

In finite volume methods, the desire to improve the order of the convergence using low order 
scheme has not attracted the attention it merits ( see the introduction of [Hj). In this context, we 
can mention the work of Martin et al. ^T], where they used defect correction method, and under 
uniform mesh to propose an implicit scheme that is second order accurate both in time and space 
and uses only first jacobian for some unsteady problems. 

The aim of this article is to develop some techniques allowing us to improve the order of the 
convergence of the finite volume solutions on arbitrary mesh conditions for second order elliptic 
problems in one and two dimensional spaces. 

We prove that, starting with a finite volume solution of order 0{h) in H^- norm, we can obtain 
finite volume solution of order 0{h?) in iJQ-norm, by using the same matrix that used to compute 
the solution . 

The heart idea used in this article is the fameous Fox's difference correction in the context of finite 
difference methods. 

The order of the convergence of the finite volume solutions, on lower schemes, depends on the 
second derivatives of the unknown solution u . 

In one dimensional space, the second derivative of u can be expanded as a combination of the 
solution itself, its first derivative and a given data. We use this idea to obtain an optimal approx- 
imation to the second derivative by using the values of the basic finite volume solution . This 
approximation allowing us to correct and to obtain a new approximation can be computed by 
the same matrix that used to compute u^, called first correction, of order 0{h^). 
Other variant to compute an optimal approximation to the second derivative of u is to use the 
fact that is satisfying the same equation that is satisfying by the solution itself but for different 
second member and boundary conditions ( this holds for some second order elliptic problems). 
This allowing to obtain an optimal approximation to the second derivative, by using always the 
same matrix that used to compute . 

We can repeat this process, successively, to obtain finite volume approximation s of orders 0{h^^^), 
where k is integer by using the same matrix of original system. 

In two dimensional space, we use the second variant that used in one dimensional space. For the 
Laplacian model, the second derivatives of the unknown solution satisfy the same equation that is 
satisfying by the solution itself, and by the same trick that used in one dimensional space, we can 
obtain a new approximation to the unknown solution of order 0{h^). 

These resultes can be extended for some Dirichlet models —Au + pu = f and to obtain corrections 
of arbitrary order we wish. Some numerical tests justifying our theoretical results are done, too. 



2. In One Dimentional case 

2.1. Basic Results and Preliminaries. The results of this article are presented in the context 
of classical functions space. We denote by C™(f2) {fl in our paper is either an interval in IR or 
rectangle in IR ) the space of continuous functions which together with their derivatives up to 
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order m inclusive are in C{Q). The norm is 

lkllm,oo,n = max (m_ax|D"u;|) 

|Q|<m O 

In all that follows the letter c stands for a generic, positive number, different at each appearance 
but 'constant' in that is independent of discretsiation parameter T,i,j 

Remark 2.1. To show that the improvement order, will be presented (in one and two dimensional 
spaces ), hold for an arbitrary admossible mesh, we try to bound each expansion with respect to 
, hi, ... (and we do so for two dimensional space ). 

Basic results given here are done in Eymard et al. "T. Let / be a given function defind on (0, 1) 
and consider the following equation 

~Uxx + auj; + /3u = /(x), X e / = (0, 1), 

(1) 



u(0) = = 0, 



where (a,/3) G ]R+ x ]R+. 

Let T be an admissible mesh in the sens of 0, i.e. given by family (-ftri)i=i,...,Ar, N G IN*, such 
that Ki = {xi_i, Xi^i) and a family (a;i)i=o,...,JV+i such that 

a;o = a^i = < xi < a;| < ... < Xi_i < Xi < x^^i < ... < x^ < XJ^_^_l = xn+i = 1, 

and 

hi = m (Ki) = Xi_^i — x^_i, for i e {1, ...,N}, 

h~ =x^- x^_i , /i+ = 1 - x^, for i e {1, ...,N}, h^ = h^^-^ = 0, 

hi^i — Xi+i — Xi,i — 0, ...,N, and size(T) — h = maxj/i^, i = 1, N}. 

The system to be solved for the three points scheme by finite volume method is 

Au'' = b, 
uo = un+1 = 0, 

where u'^ = {ui, unY, b — b^Y with A and b are defined by 



(2) 



/o\ ^ A h\ 1 / Ui+l-Ui Ui-Ui-1 \ \ , o ■ 1 AT 

(3) iAu'')^ = — + — + a{ui - + (3u„ i = 1, ...,N, 

(4) ^- = -r I f{x)dx,i = l,...,N. 

"■i JKi 

The following theorem (see |7] ) gives the order of the convergence of the finite volume solution of 
the scheme of three points (01. 

Theorem 2.1. Let f G C([0, 1]) and let u G C^([0, 1]) be the unique solution o/Ol Let t be an 
admissible mesh. Then there exists a unique solution of\^ and the error is bounded by 

(5) I E ^^"i-^l <^MM^j, 



,0<i<N 



H+4 



(6) |e,| < c/i||w2.||oo, G {!,..., AT}, 

where cq = cn+i = and Ci — u{xi) — Ui, for all i G {1, ...,N}. 



4 



B. ATFEH AND A. BRADJI 



Remark 2.2. The uniform estimation ^ yields that the order of the convergence in norm is 
at least 0{h), but numerical results shows that in general the order is 0{h^) when a — f3 — 0, this 
means that 

(7) (Y.h,^^^ <ch\ 

Before we will be able to give general formulation of an arbitrary correction, we present at 
first the first correction and after we give the second one, where additional tools will be used. 
The general formulation of corrections can be given later, by using the ideas of first and second 
correction. 



2.2. The First Correction. By integrating both sides of equation over each finite volume Ki, 
we get 

(8) —Ux{x^^i) + Ux{x^_i) + au{x^^i) — au{x^_i) + (3 / udx~ I fdx, 

^ ^ ^ ^ JKi JKi 

Then, the order of convergence of finite volume solution depends on the : approximation of the 
flux, values u{x^j^i), u{x^_i) and the integral udx. 

We shall use this idea combined with one of Fox to improve the order of the convergence of the 
basic solution on the same scheme, i.e. using the same matrix A that used to compute the basic 
solution and changing only the r.h.s of |21 

Looking now for an expansion to the error. Assuming u G C^(I) and using Taylor's formula, we 
can get 
(9) 

'■+2 '■2 

-R\^^+R\_^. Vze{f,...,7V}. 

where the following estimate holds 

(10) <c/i,+ iV3.U, Vze {0,...,7V}. 
On the other hand 

(11) = <x^) + h+ux + Si Vi e {0, TV}, 
where 

(12) <c/i+V2.|oo,Vze{0,...,7V}. 
Also 

(13) J u{x)dx:^h,u{x^) + ^Ux{xi)~-^Ux{x^^l) + Tl,'^ie{l,...,N}, 
where 

(14) \Tl\ < ch^h4u2x\\oo,yi e {1, ...,^}. 
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Substituting terms in JSJ by their expansions found in the equahties and H13|l . we can 

obtain, Vi G {1, ...,N} 

+ [- au{Xi) — au[xi^i) + (ihiu[xi) = I jax 



- ahju.^{xi) + ahl_^Ux{x^^l) - ^ {hj'^Ux{xi) - h^^Ux{xi-i)^ 
(15) - R]^! + R]^! - aSl + aSl_^ - (3TI . 

By Substituting U2x = ctu^ + [5u — / in the equation 1)15(1 we get 

U{x,+l) - u{Xi) u{Xi) - ujXj^l) \ , filU f \ f 

; + ^ au[Xi) — au[Xi^i) + pniU(Xi) = / fax 
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- ahlux{xi) + ahj_^u-j,{xi^i) - ^ [{hlYux{xi) - {h- fuj:{x^^i)) 
(16) - Rl^, + R]_,_ - aS} + aS}_^ - PT^. 

After having found an appropriate expansion of the error, we can correct the basic solution by 
approximating values and pointwise derivatives of the unknown solution u in this expansion, by 
theirs corresponding values and partial values (forward approximation ) of the basic solution u^. 

The new solution u'l = {uDfL^, called correction, obtained after these changes, will be defined 
on the same scheme, i.e. using the same matrix A that used to compute the basic solution u'^, i.e. 
Uq = ujy^j^ = and Vi G {1, N}, we have 

+ ^^^p^ + a{ul - uU) + (3h.ul ^ [ fdx 

K+1 -ht ( . ^ c\ , K - /i+ 1 



(adiUi + (3u^ - /(sj+i )) + — — ^ ' ^ (adiUi^i + - f{x^_i)j 



(17) - ah+diUi + ah+_^diUi^i - - {{hi f dim - {hrfdim-i) , 



where diUi — " 



2 

h,, 1 ■ 



2.3. The Convergence Order of the First Correction. To analyse the convergence of the 
first correction, we follow the same proof that used for proving the order of the convergence of 
the basic solution. By subtracting from lfT7|) side by side, the error e] — u] — u{xi) will be 
satisfied 



C: 6, 



where 
(19) 
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(20) Si = -| {{h+ fidiu, - u^ix,)) - - u^ix^.i))) + PTl. 

Multiplying both sides of p8|) by ej and summing from i = 1 to i = A^, we get 



^ .1 



Af 1 1 



N 



N 



■ 1 ^i+i 

1 = 1 '^2 



l2 



1=1 
N 



N 



N 



(21) 

This gives that 
(22) E 



i=l 



(el+i-eir 



hi, 1 

0<'i<A' *+2 



< 



E ^ 

^Q<i<N 



E 



AT 



— lli, 1 
0<j<A *+2 



IE '5^1 



To estimate second term in the r.h.s of H22(l . we handle each term in H20|l . Indeed, we have by 
using inequality (|14|l and reordering sum of second term of H2U|) . we can get 

N / N N N \ 

j=l \i=l i=l i=l / 

/A AT A j,i2 N 

< c I h^hf\diUi -Ua;{x^)el\ + h^h~\diUi^i - it^(xi_i)e- 1 + [^^]H^h^ef] 



i=l 
A 



A 



N 



< 



A i2 A 



1 2,1 

+ ft 

1=1 j=l 



i=l 
A 



4mE'^^(^i"»-"-(^»))']^E^^'t + mE^^"(^i"-i-"-(^-i))']^E'^»- 



.12 



But 

(23) 
and 



i=l 



AT 



\i=l 



> 0<i<A 



AT 



^A-1 



(24) 

Therefore 



, i=0 J 

- I E 'i^+il^iw* ~'«2;(a;<))^ 



, 0<j<A 



AT 



(25) \Y.Slel\<c h[ ^ +(E^1 |(E^*^' 



A 



V 



> 0<j<A 
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Inequality H25|) combined with discrete Poincare and triangular inequalities imply the following 
inequality 

^0<i<Ar '+5 / \Q<i<N 

I 0<i<Af 



i,0<i<Af 

E f^^Hi^liA + f E "'^.+^(^.^)^ 

^0<i<N J \Q<i<N 

(26) + c\hi E +fE^ 

\ \o<i<N J \i=l ' y 

To estimate the r.h.s of inequality (|26|l . we need the following estimates 

Lemma 2.1. Let = (ui) be the basic solution defined by 0), the following estimates hold 

1- I E K+^iu^~uix^+^yf \ <ch\\u2x\\ooj- 
)<i<N J 

1 

E K+^{^lU^-U^{x^^l))^ \ < c/l||'a2:r|loo,/- 
)<i<N J 

1 

3- ( E - Ux{Xi)f \ < c/l||u2a;||oo,J 

)<i<N J 

E I < ch^\\u3x\\oo,T■ 

E f'^+i^A <c/i'll"2.||oo,/ • 
)<i<N J 

Remark 2.3. Estimate 1. of lemma\^\ can be done through uniform estimate 0), where the es- 
timate holds for one dimension (see remark 2.7, page 18 in 7 ). The proof, we wish to present, 
holds also for the case of the finite volume scheme of five points in two dimensional space. 

Proof. 
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1. By triangular inequality , we have 

0<i<N 0<i<N 0<i<N 

(27) < I +^H^2.\\oo,T, 

yo<!<w J 

using the fact that cq = and X)o<i<Af ^i+i = 1 combined with Cauchy-Schwarz inequality 
to get 

yO<i<N j \^<i<N '+2 

this with H27|) imply the desired inequality 1 of lemma. 

2. Using the same thecnique, yields 

,0<i<N J \0<i<Ar 



< c/l|lM2a;||oo,/- 

3. can be obtained as done for 1. and 2. 

4. according to inequalty we have i| < cft.^, this implies the inequality 4 of the 
lemma. 

5. according to 112(1 . 5?^ is of order 0{h^) in uniform norm, which implies 5 of the lemma. 

6. according to H14|l . < ch hi, the inequality 6 of lemma will be obvious. □ 

Coming back now to the lemma ITTl since ft,^-^ — hj' = 0{h), then we have the following 0{h) 
improvement. 

Theorem 2.2. // the unknown solution w o/Q] belonging to C^{I). Then the error in the first 
correction defined by \ll\ is of order 0{h^) in the discrete Hq norm, i.e 



(28) ( E <-'^'\\-hooj, 

\o<i<N '+5 J 

where e\ = u\ — u{xi) and {u\)i are the components of the first correction defined by j_?7| ). 

2.3.1. Other Variant to Estimate the Second Derivative of Unknown Solution. For some 
cases, like the model —U2x + fiu = f , we have other possibility to approximate the second derivative 
U2x of u. Indeed, U2x satisfies the following equation 

~Vxx + /3w = f2x{x),X e / = (0, 1), 

= -/(O), 
^(1) = -/(!)■ 
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Then U2x satisfies the same equation that is satisfying by u, this allowing us to get a finite volume 
approximation to U2x, provided that u G C*{I) (see theorem I2.1|l . by using the same scheme 
that used to compute the basic solution m'', more precisely, we use the same matrix, that used 
to compute u'', to compute a finite volume approximation to U2x ■ This idea can be used also to 
compute higher order of corrections. 

2.4. Second Correction. The situation in the second correction is different to that of the first 
correction, because it is easy to pass from the derivative into its forward approximation by an 
order of convergence 0{h) (see lemma ITT|) . To get the second correction of order 0{h^), we have 
to look for approximations of first and second derivative of the unknown solution, of orders 0{h?). 
That is why, we discribe how to overcome this difhculy. 

Assuming that u G C'^{I), by similar way to that one used to compute an expansion for the 
error (21), we can get 

u{xi+i) - u{xi) u{x,) - ujx^^i) \ , r:iu t \ [ 

+ V au[xi) ~ au[xi^i) + phiu{xi) = I fax 



K+i - K I . K+i - K+iK + hf 

- ^ U2x{Xi,^l) U3x(a;,+ i) 

, K - K-l , ^ , - + K-1 I X 

+ 2 "2a;(a;i_i) + u:ix{Xi_i) 

- ahfux{xi) - a^^U2xixi) + ahf_-^Ux{xi-i) + a U2x{xi-i) 

- ^ {hf'^Ux{Xi) - h~'^Ux{Xi-i)^ - ^ {hf^U2x{Xi) + h^^U2x{Xi-l)^ 

(29) + ^hf\_^U2x{x,-l)~R^,+ ^^+R',_^^-aSj + aSU^PTl 

where 

(30) \Rl^^\ < ch',^^\\u\U^J, Vz G {0,...,7V}, 

(31) \Sf\<chf\\u\\,^^j, V*G{0,...,iV}, 

(32) \T^\<ch''h,\\u\\,^^j, V*G{l,...,iV}, 

In order to get correction of order 0{h?) taking into account the coefficients of the pointwise 
derivatives in the r.h.s of l|29|l . we have to find approximations of order 0(h?) to pointwise first 
and second derivative, and 0{h) to the pointwise third derivative in discrete X^-norm. 
Begining by the pointwise second derivative U2x{xi)j and looking for approxomation = {u2x)i^ * ^ 
{0, A^}, the idea that we want to suggest, is based on the use of Taylor's formula and values of 
the first correction. Indeed 

U2x{xi) = au{xi) + I3ux{xi) + f {xi) 

(33) = au{xi) + fU —U2x(Xt)+r^+i\+f[x^), 

where 

(34) r^+i < ch^,+ i\u3x\oo- 
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Let 5'' be the positive number 1 + ^h^_^_i, the equation 1331) becomes as 

(35) S':u,,ix,) = auix,) + (3 |^!f(£l±i)^JfM^ + /(^^) + pr,^,. 

Because of the trivial inequahty 1 < 6^ < 1 + ^, we can suggest the following approximation 

(36) (^2. )^ = —u^^^\ ^ ^ + Vz e {0, TV}. 

i ^ \ 2 / ^ 

Looking, now, for an approximation Ugj.'^ — {u'^!^)i, {0, N} to pointwisc third derivative. Because 
of 1*32; = a(3u + (a^ + /3)Ma; — /a;, it is useful to suggest the following approximation 

(37) iut% = aPul + ia' + P) |^!^k_^^ _ /^(^J, G {0, iV}. 

Remark 2.4. We can use the approximation of the second derivative, that used in the first cor- 
rection, to compute its approximation to obtain second correction. 

We shall prove now the following lemma 

Lemma 2.2. // the solution u of the equation 0) belonging to C^{I) and u^l^ '"^^x discrete 
expansions defined respectively by and j^*?! ). Then the following estimates hold 

/ N \ 3 



1- [^K+^iiu2x)^-U2Axi}f ] < c/l^ll"ll3,oo,7- 



\i=0 
' N 



2. $]/i,+ i((4'.').-"3x(x,0)M <ch\[ 



l3,oo,/- 



\i=0 



Proof. Substracting 1)36(1 from ((35|l . to get 

/ ^ r h,2. ^( ( \ i\ , P ( u{x,+ i) ~ u{Xi) U'+i-uA , 
U2x{x.) - {U2x\ = ^K^0-«,) + ^(^ K:^)^^ 

This implies , using triangular inequality with bound uniform of 5^ , that 

(N \k/N \ 3 

Yl^i+h'^^^'ix)^~U2x{xi)f\ < c[f ^/i,+ i(M(a:i) - w,i)M 



\^ ,u{xi+i) - u{xi) 
^ ^+^^ ~ 



i \2 



N 



(38) + E'^^+r' 

Using inequalities H28|l and ((34|l to get the desired estimation 1 of lemma 
By the same way, we can prove the second inequality. □ 

After having acheived optimal approximations for the pointwise second and third derivative, we 
look now for optimal approximations for {u2x{xij^i.))i, {u3x{x^j^i))i, {ux{xi)i, \fi G {0, ...,N}, we 
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have U2a; (x 1 ) = U2xixi) + hfu^x{xi) + s;, where \si\ < c/i^^||u4j;||oo, this allows us to suggest 
the folowing approximation for {u2x{Xi^i))i 

(39) uZ^'^ = (wai^), + /i+(w3;^)i, 

for pointwise third derivative, we can suggest the following approximation 

(40) =[U3xh, 

for pointwise first derivative, we can use the trick that used for pointwise second derivative 

/X u{Xi+l) — u(Xi) ^l+i , , 

(41) uix,) = ^ ^ - -^U2xix^) + U, 

where \ti\ < c/i^ i ||m32;|oo, and an approximation will be suggested as follows 



(42) (-.^'^). = %^-%i(4f). 



We would now prove the following lemma 

Lemma 2.3. If the solution u of equation 1 belonging to C^{I). Then the approximations (u2^^'^)^q, 

('"3x^'^)ilo '^"■'^ ii^x''^ )i)iLa defined respectively by the expansions f39\) . and satisfying 

the following estimates 

/ N \ i 



/ ^ 1 \ ^ 

2. \^^hi^i{ulp'^ -U3x{x,+ i)f\ <ch\\u\\4^^^j. 
3- (E/i,+i((«L''').-u.(a;.))M <ch' 



'^'|w||3„,/- 



\ j=0 

Proof. 



1. We proceed as done in the proof of lemma Using triangular inequality and equality 
(|^ . to obtain 



N \ 2 / N 

i+i,2 



'^'^'^+^iuZ^' ~ U2xiXi^l))^ ^ < [J2h,^l{iu2x)^-U2xix^))^j 

N 

h,+ ih+ {{u'^-^)^ - U3x{xi)f 

^ AT 



^1=0 / \i=0 

/ N 



Ki=0 



Using lemma 12.21 and the uniform bound of s,; to obtain the desired inequality 1 of the 
lemma 12. 31 

2. and 3. of the lemma can be handled by the same way as done for the first estimation. □ 



12 



B. ATFEH AND A. BRADJI 



Now we are able to define the second correction = (wf)^^^, wliere Uq — uj^_^_^ — and for 
all i e {1, N}, we have 



^i+l 



+ — au^_-^ + (3hiU^ = / fdx 

, K - K-1 »-i2 , -Kf^t-l+'^t-l^ i-i,2 

1, + / /i.2\ "i / h.2\ I ,+ / ft.2\ ( "i-l / h,2N 

- ^{hf{<% - hfiu'i^.,} - ^{ht'{u^\ + V'(4.').-i} 

(43) + ^hfK.^{u\%., 

To analyse the error of the convergence, we proceed as done for the basic solution and the first 
correction. Indeed, let ef = uf — u{xi) and substracting equality from to get 

(44) _f|t^ + + aei - aeU + Phef = - + S'f- 
where 

7» = ^ -W2x(a;,+ i)) -U3a;(a;i+i)) 

(45) - ah+i{u'^-% - u^ix,)) - c.^{{u\^\ - u^^^Xi)) + R\^. + 



= -|{/if ((«^'^).-^..(a:,))-V'((u^'^).-i-z..(x._i))} 

(46) - ^/^-\,_l((«t2)._l-^i2.(a;.-l)) + /?I^2 

Using the proof of the convergence of the basic solution and the first correction together with 
triangular and discrete Poincare inequalities combined with lemma and expansions (|30|l . H31I) 
and l|32|) to get the following 0{h^) improvement. 

Theorem 2.3. // the unknown solution of 0) belonging to C^(I). Then the error in the second 
correction defined by \4-S^ is of order 0{h^) in the discrete Hq norm, i.e. 

(47) I E <c/.^lkllw 

\0<i<N '+3 J 

where ef = uf — u{xi) and (uf )i are the components of the second correction defined by 
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2.5. Corrections of Higher Order. In this section, we give the general formulation of an arbi- 
trary correction. The pointwise derivatives will be approximated in the light of ones of the first 
and second correction. The proof of the order of the convergence is the same one that done for the 
first and second correction. 
For each integer k > 1 

u{xi+i) - u{xi) u{xi) ~ uix^^i) N , / \ , i / ^ 
+ = -u^(Xi+i) +u^{x^_i) + a^u„,^(x.^^i) 

'''^2 '^~2 'in— 2 

m=2 

and 

(49) <c/i^+i|klU+2,oo,7 
We have also 

(50) u{Xi+l) = U{X^) + E ' I Umx{Xi) + S'f , 

Tfl. 

m— 1 

where 

(51) \S^\<ch^.^'\\u\\,+,,^j 
For the fifth term in 0, we have 

(52) / u{x)dx = h,u{x^)+Y^ 7I7--ni"™^(^')~Xl H ^7p^W(m+J>(a;^_l)+T;^ 
where 

(53) \T^\<ch.h^^'\\u\\^+,^^j 
Substituting terms of © by theirs expansions H50|l . (|51|l and (|53(l . we obtain 

T \ r + au{xi) - au{Xi-i) + phiU[x^) = / /dx 

fc+i fe+i 

m— 2 m— 2 
- OL y T—U^,j:{Xi)+a } j Umx[Xi-l) 

^-^ m! ^-^ ml 

m—l m—1 



^-^ m + ll! ^-^ (m + 1)! ^-^ v. 

m=l ^ ' m=l ^ ' j=0 ■' 

(54) - 1 + i - "^^'^ + "^t 1 - 



After having found an expansion approximating the equation (jSJ , we need now the following useful 
lemma 
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Lemma 2.4. Each kth derivative of the solution u of the problem^ can be expanded as a linear 
combination of the solution itself, its derivative and the derivatives of the given function f up to 
and including ik — 2)th derivative, i.e, there exist reals {cKjjj^o {^i'^2} such that 

k-2 

(55) Ukx = X! /j^ +a'lu + a!^u^. 

j=o 

Proof. We can prove this lemma by induction on the integer k. □ 
Assuming, now, that we have obtained the (k — l)th correction u^ -^ = {u^~^)i, i.e. approximation 
of order 0{h^). According to equahty H54|l . to obtain correction of order 0{h^'^^), we have to find 
approximations for the pointwise derivative up and including fc + 1 order of the solution u. The 
idea which we will present is similar to that one presented to compute second correction. 

At first, we look for optimal approximations to U2a;(a;i),..., W(fc+i)2;(a;i)- To do so, we use the 
previous correction, i.e. (fc— l)th correction, and the optimal approximations to U2x{xi)^...^Ukx{xi) 
used to define this correction. That is why, we define the kth correction by induction, we assume 
that, we have obtained {k — l)th correction of order 0{h^) and we have found optimal approxi- 
mations (according to theirs coefficients in (|54|l (u2i'^^^)i,---i(Wfc^'^~^)i for U2x{xi)^...^ Ukx{xi), i.e. 
theirs orders of convergence are 0{h''~^),...,0{h). 

Because the coefficient of U(k+i)x{xi) in (|54|l is of order 0{h'^), it suffices to approximate it by 
in order 0{h). This, can be done easily through lemma i.e. an approximation defined by 

fe-i 

(56) = ^af V,.(^.) + «^'"r' +«2+'5"'"\V^ e {0,...,iV}. 

0=0 

We can use, also, in H56(l instead of the the (fc — 1) correction, the basic solution u^. 
For any integer (j such that 2 < (3 < k, we look to find approximation u^^ of order 0{h^^'^~^) to 
pointwise derivative {upx{xi))i of order (3, bacause the coefficients of such derivative in H54|l are of 
order (3 — 1. We have through lemma [T^ 

13-2 

uisxixi) = a'^fjxjxi) + a1u{xi) + a^Uxjuj) 

3=0 

0-1 

= fjxjxj) + gf M(xi) 



(57) 



3=0 

( , . , . k-P+2 



J=2 J 



where 

(58) kfl <c;if+J-''|uU+3-/3,oo 

An obvious approximation for pointwise derivative of order (3 can be given as 

P-2 - l/fe-l '=-/3+2 /, fc>fc-lv 

(59) i:^^). = E «f /.^(-o + ^f-r^ + 4 ^-'^ — E 

3=0 V j=2 



We would prove the following lemma 
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Lemma 2.5. If the solution u of the equation^ belonging to C'"'^^(/). Then the approximations 
ujj^ , where 2 < /3 < + defined by the expansions \5f)]) and i|5.9|] satisfying the following estimate 

(60) (^K+^S^u';'),-U0,{x,)f^ < ch^-^+^\\u\\^+,^^j. 

Proof . We can prove this by induction. □ 
After having found optimal approximations to fundamental pointwise derivatives, we derive now 
optimal approximations for (u^;j;(xj_|_^))^g,2 < (3 < k + 1 and {ux{xi))f^Q. 
We have 

(61) ui3x{x,+ i_)= ^^U(/3+J>(a;^) + s^V^G {0,...,7V} 

3=0 ^' 

where 

(62) l4l<c(/»+)'=-'^+'k(fe+2).U 
We can suggest the following approximation 

(63) (</'').= E %^(«J>'+,.)J.,Vze{0,...,iV}. 

3=0 •'■ 

For the pointwise first derivative, we can do 

(64) Ux{Xi) = — —Ujx{xi) + ti, 

'^^+h j=2 ^■ 

where 

(65) <c/i^+i|w(fc+i)x|oo,V*e{0,...,iV} 
this allowing us to consider the following approximation 



fc-i _ fc-i k \ 
*+2 j=i ■' 



(66) {^^% = - E Vz e {0, ...,iV} 



We need the following useful lemma 

Lemma 2.6. If the solution u of the equation ^}) belonging to C*''+^(J). Then the approximations 
("/3a:^ where 2 < /? < fc + 1, defined and (u^''"')i defined respectively by the expansions ^6y\) and 
satisfying the following estimation 

1 

/ Af \ 2 



AT \ ^ 



\i=0 

/ N 



\i=0 
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Proof. The proof can be done as done for proving lemma □ 
Now we are able to define the kth correction wjj — where u§ = = and for 

i e {1, N}, we have 



+ — + aUi - au'l_^ + phiUi = / fdx 



fe+l k+l k k 

Z / Z / ^ / jjil Z / ,^1 

m— 2 7n—2 m—1 m—1 

Consedering the following expansions 

(68) 7f = - E )^ - - « E H^^^""'^^ " ""^^^"^^^^ + ^^+^ + 

rn—2 m—1 



(69) 



Let = — M(xi) be the error in the kth correction, thus 

(70) _fitLZ^ + + ae'l - ae^ + (iKe\ = 7^ - 7^-1 + 

Using the proof of the convergence of the basic solution, first and second corrections and lemma 
Eini together with (gHl jEHIi, (ED and to get the theorem 

Theorem 2.4. If the unknown solution of 0) belonging to C^'^'^{I). Then the error in the kth 
correction, defined by j70| ) is, of order 0{h^^^) in the discrete Hq norm, i.e. 



u 



fc+2,00,/ 



(71) ( J2 ^^^^^1 <ch''+' 

\o<i<N "-i+h J 

where e^ = — u{xi) and (w*^); are the components of the kth correction defined by | |6'7| ). 

Remark 2.5. As you have seen, we can generalize the results obtained for general equation y" 
f{x,y,y'), where f is a smooth function. 
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(72) 



3. In Two Dimension Space 

3.1. Basic Results. Considering the second order elliptic problem, with homogeneous boundary 
conditions 

-Alt = f,onn = (0,1)2 

where F = d^l is the boundary of and assuming that the solution u is belonging to C'^{Cl), and 
the second member / e C(f2). 

Let T = {Kij)Ki<M;i<j<N be an admissible mesh of Q in the sense of \7\, that is satisfying the 
following assumption 

Assumption. Let M, N E IN*, {hi)f£^, (^j)jLi '^'^e positive numbers and such that 

M N 

j2hi = Y,kj = 1, 

z=l j=l 

and let Hq = Hm+i = h) = ^'w+i = 0. We define: 

XI = 0, for i G {1, M} : a;,+ i = i + hi, 

yi = 0, for j e {1, ...,N} : y^+i = yj_i + kj, 
(So that Xm+i/2 = yN+i/2 = Ij; and 

Kxi = [xi_i,Xi^i], Kyj = [yj-_i,yj+i], K^j = Kxi x Ky,. 

Let (xi)^^^ and {yj)^^^ he points such that 

Xi_i < Xi < Xi_i, for i — 1, M; xo — 0, xm+i = 1, 



y,_i < yj < y.-i, for j = 1, ...,N; yo = 0, yN+i 



1, 



and let Xij — (xi,yj) for i ^ 1, AI and j — \, N . Set 

h^ ^ Xi ~ Xi^i, h'l = x^^i - Xi, for i = 1, M, ft-j^i = Xi+i - Xi, for i = 0, ...,M, 

kJ ^ Vj - yj-i , k+ = y^+i - yj, /or j = 1, ...,N, k^^i = y^+i - y^, /or j = 0, ...,N. 

Assuming that Hq — /ij/^x ~ ^o^ ~ ^iv+i ~ and considering the mesh size h = max{{hi,i 
1,...,M), (fc, = 1,...,7V)}. 

Definition 1. Let X{t) be the set of functions from D, to TR piecewise constant over each Kij. 
For w € X{t), we define the discretes H^-norm and L^-norm respectively 



(73) 



I 

. i = 0,M 

\j=0,N 



1/2 

/ 



Tnrr -J^ 2 



(74) 



\\w\\l2 



. .= 0,M 



\ 



\ 



Let w = {wij)o<i<M+i.Q<j<N+i and be the the following discrete opertor 



(75) (A^w), 



Wi+l, 



Wi. 



Wi. 



Wi-l,j 



+ l ~ Wi 
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for a continuous function g £ C(f2), we introduce a similar definition: 
(76) 

(A g),, - -k,[ — — )-h,{ — - ^— ). 

/■-r 2 2 ^^^2 ^^2 

To simlify the notations, denotes ^ i^-rn Sij denotes ^ i=o~M" 

3=TN j=0,N 

Integrating the equation (|72(l over each finite volume Kij , to get 

(77) -[ (..fe,.ri-»,(x._,„,))%-/ 
Taking the first term in left hand side (l.h.s) of (|77|l 

where 

1 

(78) i j= - / {y -yj)u^y{x,+ i,yj)dy, 

J y . 1 

and is some point lies between yj and y. 
Using, again Taylor's formula, yields 

u{xi+i,yj) - u{xi,yj) 
Thus the following estimates hold 

(80) < cA:||Mxj,|oo,f2> and i jl < ch^+^\u2x\oc>,n- 

By the same way, we can get 

(uy{x.y^^,)~uy{x.y^_,))dx = -K ^(^^' ^.-+0 ' "(^-- _ <^^.V,) - <^^.V,-^) 

(81) -hi - - ^1,3+^ + 

where 

(82) I'S'.j+il < c/i-|uxy|oo,n> and \R^J+^ \ < ck^+^\u2y\oo,n- 
Therefore, the equation 172|) becomes after integration as follows 

(A*^^u)y = / j-i?,_i^j)+/ii(i?,j+i-i?ij_i)+S',+ i j-S',_i j + 

J Kij 

The basic finite volume solution — (uy )i=o....,A/+i j=o,...,a'+i is defined by 

(83) uoj = UM+i.j = Uio = Ut^N+i = 0, 
and for (i, j) G {1, M} x {1, N}, we have 

(84) {A^u% = f fdxdy 

The existence, uniqueness of the solution , the analysis of the order of the convergence can be 
justified as done for ID case (see [7]). More precisely, we have the following theorem 
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Theorem 3.1. fjTj^ // the solution u of the equation J7i^ belonging to C^(f2) and f G C^{C!,). 
Then the approximate solution — (uij) defined by the boundary condition l\8S\) and the discrete 
equation satisfies the following estimates 

(85) l|e||i,r <c/i||w||2,^,o, 

(86) Ml- <ch\\u\\^^^^^, 

(87) (^k,h,+ .e^^}j <ch\\uh^^^n- 

where Cij = u{xi,yj) — uij for {i,j) £ {1, ...,M} x {1, ...,N} and vanishes elsewhere. 
3.2. First Correction. In this section, we assume more regularity for the solution u, i.e. u G 
Looking, again, at the equation l|77|l . to simplify the notation, let 

9t{y) = u^{xi+i,y) - u^{x.^_i,y). 

Using Taylor's formula, we get: 

(88) - / 9i{y)dy = -kjg,{yj) {gi)y{yj) - Tij, 

where 

2y I oo,Kj 
< C.k)h,\u2x,2y\oo,n- 

Using again, Taylor's formula, yields 

{k+^ - kj\gi)y{yj) = {k+^ - k-\u,xy{x.^^i,y.j) ~ Uxy{x^_i,yj)) 

(90) 2 _2 

= A:+ hiU2x,y{xi,yj) ~ kj hiU2x.y{xi,yj-i) + 11}^ 



k/)S}^ - {kf - kr')Sl 



where 



(91) K\<chf\mx, 



(92) \Sl\<ch- |«3.,,loo,a 

(93) \Ulj\ <ckT^h,k^_l\u2x,2y\oo,n 
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By substituting gi by its value in l|88|) and by using equality (|90l) . we get 

{u^ix,+ i,y) - u^ix,_i,y)j dy = -kj (^u^{x,+ i,yj) - u^{x^^i,yj)j 



(94) -hi^^U2x,yiXi,yj) + h,^-U2x,yiXt, yj^l) 

but, in the other hand 

u{xi+i,yj) - u{xi,yj} K+i - '4 , ^ i 7,2 

(95) =w:r(a;i+i,yj) + ^ u^^ix^^^^y^) + R^^^_^^, 

where 

Combining and (|^ yields that 

^„ _ f ^^+i,yj) -uix,,yj) u{x,,yj) -u{xi-i,yj) 



(ux{Xi^i,y) - Ux{xi^i,y)j dy 



u 1 - - \ ^j+i ^i-i / 

tfj-i \ '+2 ' 2 / 

_^ , / K+1 - K , . K ~ I 

^ \ 2 "2x(a;i+i,yj) ^ M2a;(2;,_ 1 , %) 

k ■ k ■ n 

- hi^^U2x,y{Xi,yj) + hi^^U2x,y{Xi,yj^l) +kjR^^i^. 

k+^-k^'^ 

(97) - k,R^^^- ' 2 ' ^+ ' 2 ' 

and by the same way, we can find similar expansion for the second term in the l.h.s of H77|) . 

u{xi,yj+i) - u{xi,yj) u{xi,yj) - u(xi,yj-i)\ 



^ (uy{x,yj^i) ~Uy{x,y^_i)^ dy = -hi 



n, ■ - fc^ k~ — k'^ 

hi ( -^ii^ -U2y{xi,y.j^i) - -^—Y^U2y{xi,y^_i) 



- kj-'^Ux,2y{Xi,yj) + kj-!^Ux.2y{Xi-l,yj) + 

(98) - h^K^^. - - + 2 ^^' "-^ - ^ - 

where 

(99) ^c^j+il^a^'loo.n 

(100) =^cfc+Vx,3y|oo,0 

(101) =^cfc~V.,3yLn 
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(102) \H,,\<chlk,\u2.,2vU^ 



(103) \U!^j\<ch- kjhi^l\u2x,2y\oc,n 

Equalities and combined with (jZH) yields 

h7. , ~ht . . . hV - h1 

fdxdy — k, 

Ik, 



(A*^^u)y = f fdxdy- kj ^ U2x (x,+ 1 , ) + % — — ^^U2a,{xi_i,yj) 



-hi^^ — -U2y{x^,yj^i) + hi^ — 7r~^U2y{xi,yj^i) 



^-2 ^_^2 ^_2 

+ ht-^U2x,yiXt,yj) - hi^^U2x,yiXi,yj-l) + kj-^Ux,2y{xi,yj) - kj-^Ux.,2yiXi-l,yj) 

7+2 _ , -2 ,^2 ^ , _2 

(104) . + . + ^^'^^^' 4 - + i?.. + ^ + c4- 

Let Ti = [0, 1] X {0}, = {1} x [0, 1], = [0, 1] x {1}, r4 = {0} x [0, 1] be the partial boundaries. 
The second derivative U2x of u is the solution of the problem 

-Ai; -/2,, on r!= (0,1)2 

^^Iti = w|r;, = 0, 
^^Ir2 = -/|r2, 



(105) 

Let v'^ = (vij) the finite volume approximation of U2x, then 



ViO = Vi^N+l = 0, 

(106) { voj=-f{0,y,), 

VM+1,3 = -/(1,2/i) 

and for (i, j) e {1, M} x {1, N}, we have 

(107) {A-'v^, = / /2.d2:dy. 

Remark 3.1. We remark that, we have used the same matrix, that used to compute , to compute 
an optimal approximation for U2x- This implies, in turn, that we have the same order of the 
convergence as for the solution u provided that u belonging to C'^{Vl), i.e. theorem \S.l\ holds also 
for U2x instead of u and instead of u^. 

Therefore 

Lemma 3.1. Let u he the solution of the equation jyij| ). If u £ C^(f2). Then the approximation 
defined by the boundary conditions U06]) and the equation \10T[ ) satifies the following estimates 



(108) \^hikj{u2x{Xi^i,yj)-v^jf\ <ch\\u\\ 



4,oo,n 
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(109) 



'^h^kJ^l{u2x,y{^^,yJ) 



< c/i||w||4,oo,n 



(110) 



< ch\\u\ 



4.oo.f2 



Proof 

1-Using triangular inequality combined with estimate H86|l of theorem 13. II yield 

'^h,kj{u2xix,+ i,yj) - V^j)"^ < '^h^kj{u2x{Xi+l,yj)~U2x{x^,yjf 



'^hikj{u2x{xi,yj) -VijY 



< c/l(||u3:^||oo + ||u||4,oo,n) 



< chllul 



4.oo,n- 



2-Estimates ((TU^ and ((TTn|l can be proven by the same way, i.e. we use triangular inequality and 
theorem 13. II P 

Let = (wij) be a discrete function, to simplify the natation, we define the following discrete 
operators 



and 9f^,,^ ^'-^+^~^'^ 



Now, we are able to define a new approximation Ui — {uj^), called correction, of order 0(h'^) 
(as we will see) defined by the boundary conditions H83(l and the following discrete equation for 
(t,j)e{l,...,M}x{l,...,7V} 



(111) 

where 

(112) 



{/Vu'Dij = / fdxdy+ (7i+i,i - 7i-ij) + - 7j,i-i) + 



7.+ i , = -k: 



(113) 



^7 _ ^ + 

7ij+i = -hz^^ — ^(-/(a;,,yj+i) - Wy), 



(114) 



Sij = hi^'—dfvtj - hi^^—dfv^j-i + kj^—{-fx{xi,yj) - dfv^j) 



-kj^{-fx{x.,.i,yj) - d^v.^ij). 



Remark 3.2. As you can see that, we have used only tha approximation of U2x cind the fact that 
U2y = — / — U2x to approximate all the higher pointwise derivatives in the r.h.s of 
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To analyse the error, let ej^ — uj^ ~ u(xi, yj) be the error in the first correction, by substracting 
(|104() from we get 



(115) (^^e )y = 7,+ 1 - + - + (5, 



where 



(116) 7i+i J = -kj-^^ — -{vij ~ U2x(a;,+ i , %■)) + kjR^^i -, 



— 

(117) = — '-{u2xixr,yj+i) - n,j) + hiR^j^ 



'3 + 1 "'j 
2 

2 , _2 , _2 



fc+ kj fc: 

+kj-^{u3xixi,yj) - d'fvij) - kj-^{u3xixi^i,yj) - dfvi^ij) 

; +2 J -2 ,+2 , _2 ,2 , _2 , i 2 , _2 

ol , -fej o2 ^+ -fe. .1 I -fe. .2 

2 y 2 2 2 -' 

(118) -T,, - H,, - C/i - C/2 

3.3. Convergence Order of the First Correction. Multiplying both sides of H115|l by ejj and 

summing over e {1,...,A/} x {l,...,iV}, to get 

''"'"2 l,j 2 ■^"f'2 ij' J 2 



(119) = ^ 7»+i j4 - ^ 7^-1^4 + - II + H 

*,i *,i i,3 i,j 

Reordering equation pi9|l and using the fact that ej^ vanishs on the boundary mesh points, to get 



1 . 



(120) We'Wl, = - 4) - Il7.,,+i(et-+i - 4) + 11 ^^^4 






We should now estimate each term in the r.h.s of (|120|) . The first and the second term can be 
handled by the same way. Hence, it is suffices to estimate the first and the last ones. Using 
Cauchy-Schwarz inequality yields that 

lE7.+ ..(e,V.-ei.)l < (E%^)^(Efc.%^)^ 
(121) < (E^^T 

/Co 

Using triangular inequality, in order to get 

(E — X — ^ - E^ — — K+i{v^,'^u2x{x,+i,yj)r\ + \}_^ — 
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Using estimate (|1U8II of lemma 13.11 (combined with estimate (|87ll of theorem 13.11 ) and inequahty 
(EH), to get 



'■,3 



3,oo,Sl/ 



(122) < ch^\\u\ 



4,oo,f2 



Comming back now to the last term in the r.h.s of (|120|l . using triangular inequality yields that 

1 1 1 k 

< -^2x,yixt,yj)\\elj\+Y^ h,^\dlv,^j-i - U2x,y{xi,yj^i)\\elj\ 



ij i,j i,j 

«,i «j' ij' 

(123) + EV^-ii4i+EVi^ii4i 



Begininig by the first term in the r.h.s of H123|l . using the Cauchy-Schwars inequality and the 
estimate pU9|l of lemma 13.11 



(124) < ch?\\u\U^^jAY^h,k^+i{e\^ 



^2 



This implies that 



,1 



(125) E '^^^l^i^' -"2x.y(x,,%)||e,y-| < c/i2||e^||i^r||u||4,oo,n 
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and by the same way, we can handle the second term in the r.h.s of H123|l . Indeed 



a k 



_2 



hi^^\dfvi^j_i - U2x,yixi,yj-i)\\elj\ < hik- \dlvi^j-i - U2x,y{xi,yj-i)\\e]j\ 



( y 



L2 



Kkj^i{dfvij-U2x,v{xi,y])y 

J=OJV 



This with estimate (|1U9II imply that 

_2 



(126) 



a k 



h^^\^lv,^J-l-U2x,y{x^,yj-l)\\e\^\ < c/i^ ||el 1,^ ||w|| 



4,00,0 



Taking, now, a look at the other kind of terms. Using estimate (j91(l to get 



\Slj\\elj\ < c||u||4^^_oXl 

id 



i 



< c/i^l|w||4,oo,n 



2J 



Hence 
(127) 



By the same way, we can find the same estimate for terms corresponding to Sf^ , L\j , Lf^ , Tij , U}^ 
and Ufj . 

Inequalities ((n5|l - lfT77|) yield that 

(128) ^c'^'ll^ll4,oo,nl|eii,r. 



Combining equality H120() with inequalities H122() and (|128|) yields 

Theorem 3.2. // the solution u of the equation belonging to C^(j7) . Lei w'' 6e the basic 
finite volume solution oj boundary condition \8!^) and discrete equation \84\j . Then the finite volume 
approximation u\ = {u\j) defined by the boundary condition \8S6 and the discrete equation 
satisfies the following 0{h) improvement in H^-norm 

(129) \\e%,r<ch^uh^n, 
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(130) 



Mi? 



< ch^ 



l4,oo,n! 



(131) lle'llL^ <c/i'lk||4,oo,n, 

where e\^ — u{xi, yj) — uj^ for G {1, M} x {1, N} and vanishes elsewhere. 

Remark 3.3. Numerical results schows that the order of the convergence in norm even of the 
basic solution is 0{h^), i.e. the order in and 1^86]) is 0{h^) hut the coeffecient of the order 
in f J.WJI) and U31\l are smaller than to those ones of iS6]) and J^7| ). 

3.4. The Second Correction and Higher Order of Corrections. In this subsection, we give 
an idea allowing us to construct second correction, i.e. the order of the convergence is 0{h^), this 
result can be extended to construct an arbitrary correction we wish. To compute the second 
correction, we use the first correction to estimate some pointwise derivatives of the solution u 
and the fact the equation l(7^ satisfying by u. It suffices to remark that, provided that at least 

u e C^{n) 

^ fdxdy = -kj(^Ux{xij^i,yj)~Ux{Xi_i,yj)^ -h^(^Uy{xi,yj+i)-Uy{xi,y^_i)^ 

u J ~ J I \ 1 K ^ I \ 

-hi ^ U2x.y[Xi,y-j) - kj Ux,2y\Xi,yj) 



hf~h7\ hf-hj^ 



2 -kj^^^fxyix,. %) + Omh,) - 0{htk,) + Oi^ht) - Oik^hf). 

To approximate the pointwise derivative U2x,yixi,yj) , we compute the first correction to the 
unknown solution U2x (because it is satifying the same equation that is satisfying by u, i.e. 11051 
this first correction is of order 0(/i^) in Hq, this means that U2x,y{xi,yj) can be approximated by 
an 0{h^). By similar way, we can approximate U2x,y{xi,yj). 

Comming back to Ux{Xi_f_i, yj) — Ux{xi_i, yj), the derivatives will appeared here in the approx- 
imation of Ux{x^j^i,yj) and Ux{x^_i,yj) are similars to those obtained in one dimensional space 
for the second correction, and consequently, we can use the approximations just obtained to U2x 
to approximate such derivatives. 

3.5. Some Extensions of the Results. So far, we have considered the Laplace model, where 
the second derivative of the solution u are also solutions of the same equation. In this section, 
we attempt to extend results obtained, to some second order elliptic problems, where the second 
derivatives of the solution u are also solutions but for second member depends on the solution u 
itself, its derivatives and a given function. 

The idea will be used is to approximate these terms by theirs ones corresponding in the finite 
volume solution, i.e. u and derivatives of u will be replaced by u'* and divided difference of 
respectively. 

Let us consider the following model 

-Au + pu = /,on n ^ (0,1)2 
u|r -0. 

where p is a given function and p > 0. 

We use the same scheme that used for Laplace model. As done above, we look, at first, for the 



(132) 
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finite volume solution ( basic solution), after, we look for a convenient expansion for the error, 
where we try to approximate the derivatives of the unknown solution u by using the basic solution 
u^. 

3.6. The Finite Volume Approximation ( Basic Solution). We use the same notations that 
used in the second section, therefore, for f,pG €^{^1) 

{A*-^u)ij+hikjPiju{x^,yj) = / fdxdy + kj{R.^^i,j-R^_i,j) + h,{R^.j^i-R.^.j_i) 

J Kij 

(133) + S^^i j ~ Si^i j + S'ij+i - S,j_i - Nij, 

where ^p{xi,yj), Ri+ij, Rij+i, are defined as in (jZHjl , (CHI) , JHOll, JHOI), l|H2)l and (jH21l- 

Nij is defined by 

(134) N,, ^ ({x- x,)^(a,,) + (y - y,)^(%)) dxdy 

where is a some point in Kij . Then the following estimation holds 

(135) |iVy| < c/i,fcj(/i, + A:j)hlli,oo,n- 

The basic solution — (itij)i=o,....M+i,j=o,...,JV+i7 which will approximate the solution u of the 
equation (|132|l is defined by 

(136) uoj = UAZ+ij = UiQ = Ui,jv+i = 0, 
and for [i,]) G {1, M} x {1, N}, we have 

(137) {A''u''-)ij + hikjPijUij ^ / fdxdy. 

The existence and uniqueness can be done by using the same techniques in ID (see [J])- 
Using techniques that used in the second section yields (given in 

Theorem 3.3. If the solution u of the equation belonging to C'^(Vl),the coefficient p belonging 

to C^(f2) and f G C(f2). Then the approximate solution = i^ij) defined by the boundary 
condition Ulclb'j) and the discrete equation jl!^'/^ , satisfies the following estimates 

(138) \\e\\i,r < ch\\u\L^n, 



(139) lzZ^J^*+^4| ^ c/i||u||2,oo,n: 



(140) \\e\\L^<ch\\u\\,,^,a, 

where Cij — u{xi,yj) — Uij for {i,j) £ {1, ...,M} x {1, ...,iV} and vanishes elsewhere. 
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3.7. The First Correction. We proceed as in the third section, we begin by finding an expansion 
of the error. Begining by J^^ pudxdy 
We have 

(141) J pudxdy = hikjPiju{xi,yj) + kj— — ^ ^' + — 2 ^ '^g'^'^ ("'J^ ^^i' 
where = {xi,yj) and 

(142) \Nl,\<cMk,h:'\\u\\^^^^^. 

Combining equahties (|104|l and H141|l yields 
(143) 

( A*-''^u)ij+hikjPiju{xi,yj) = / fdxdy 

JKij 

-kj '^'+\ U2xix^+i,yj) + fc /'' J''~' M2a;(a^^^i,yj) " ^"'^ 2 U2y{xi, y ^j^l) + h,-i-^^U2yix., 

^+2 ^-2 

+ hi^U2x,y{Xi,yj) - hi^U2x,v{Xi,yj-l) + kj^Ux^2y{xi,yj) - kj-^Ux,2y{Xi-l,yj) 

-fcj ' 2 ' {Px{xi,yj)u{xi,yj) +PijUx{xi,yj)) - h, ' ^ ' {py{xi,yj)u{xi,y.j) + pijUy{xi,yj)) 

We look, now, for an approximation to the second derivative U2x of u by using the same matrix 
that used to compute the basic solution u^. Remarking that U2x is the solution of the following 
equation 

(144) -Aw +pv = f2x - P2xU - 2pxUx 

with the boundary conditions 



v\ 


Ir, 


= ^^Ir3 = 0. 


v\ 


Ir. 


= -/|r2, 


v\ 


Ir. 


= -/Ir.. 



(145) 

Hence an approximation = (vij) to U2x can be defined as 

{ViO = Vi.N+l = 0, 
vo, = -f{0,y,), 
VM+1,3 = -/(1,2/j) 

and for (i, j) e {1, M} x {1, N}, we have 
{A''v'^)ij + hikjPijVij / f2xdxdy 

JKij 

(147) - hikjP2x{xi,yj)uij - 2hl kjPx{xi,yj)dluij - 2h:[ kjPx{xi,yj)dfui^ij 

where the components of the finite volume solution u'* defined by 1136|l - (|137|l . 

To analyse the convergence of the finite volume approximation v^, let v — U2x and using equality 
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l(T^ to get 

(A*'^U2x)y + hikjv{xi,yj) = j {f - p2xU - 2p^u^)dxdy 



(148) + S,+ l^jiu2x) - S,_l^j{u2x) + S,j+l{u2oc) - (U2i.) - N,^{u2x), 

where R^_^_i .j{u2x), -Rj.j+i ("22;), Si^i j{u2x), S^ J_^_l{u2x), N,j{u2x) are the same previous expan- 
sions by substituting each u by U2x- 
This imphes that 



{A*'''u2x)ij + hikjv{x^,yj) ^ / fdxdy 

- hikjP2xixi,yj)u{xi,yj) - 2kjhfpxixi,yj)ux{xi,yj) 

- 2kjh:^px{xi,yj)uxix^^i,yj) + 0{hik.jh) + 0{h^kjh^_i) 

+ kjR,+ ijiu2x) - kjR,_l j{u2x) + h,R^ j^l{u2x) - 1 (U2x) 

(149) + 1 j(u2x) - S.^_l.jiu2x) + S^j^l{u2x) - S,.j_l{u2x) ~ N,j {u2x) 

Let — (rij)ij = {v{xi,yj) — Vij)ij be the error in the approxomation (|147|l . Substracting 1147|) 
from ifTI^ . we get for e {1, M} x {1, ...,7V} 

(150) (A^r)y + hikjp.jr.j = a^^i ,^ - 1 + a^^^.+ 1 - a.^^i + 5,j 
where 

(151) "i+i J = fcj^j+i j(w2:e) + S,+ l^j{u2x ), 

(152) a^j+i = h,R,j^i{u2x) + Si^j+i{u2x), 

{kj = -hikjP2x{xi,yj){u{xi,yj) - utj) - 2kjhfpx{xi,yj){uxixi,yj) - dfuij) 
-2kjh- px{xi,yj){ux{xi^i,yj) - dfui^ij) + 0{hikjh) + 0{h^ fcj7ij_i). 
Multiplying both sides of 11501) by and summing over £ {1, M} x {1, N}, to get 

(154) llr'^lli,^ + h,kjp,jrl = - ^ a^^ij{n+i,j - r,y) - ^ a^+i ^(rjj+i - ) + ^1 



Using the tricks those used to bound the error in the first correction (subsection 3.3) to obtain the 
following optimal approximation to U2x 

Lemma 3.2. Let u he the solution of the equation ^I'JSJjl . If u & C'*(f2) and p S C^(f2). Then 
the approximation defined by the boundary conditions il4(>V o,nd the equation J j/^7| ) satifies the 
following estimates 

(155) \Yhikj{u2x{xi+i,yj)-v,jf\ < c/i||w||4 



(156) 'Yh^+ikj{u2xix,+ i,yj) - v^jf < ch\[ 



l4,oo,n 
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(157) '^kj+ihi{u2a,{x,^i,yj) -v^jf < ch\[ 



l4,oo,0 



(158) I J2h,k^_,.{u,,,y{x.,y,) - !^M±L_^)2 I < ch\\u\U^^^^ 

■i J+ 2 



(159) |^fe,^.fc,-(^3.(a;.,y,)- "'+;^'^^/'-' f I <ch\\u\ 



4,oo.f2 



After having acheived an optimal approximation for U2x , we have to find an improvement of the 
basic solution u'*, i.e. correction of order 0{h^). Looking again at the equality (| 143(1 and rewrite 

(A*'^^)^^ + hikjPiju{xi,yj) ^ [ fdxdy-kj — -U2x{xi+i,yj) + kj— — —^U2xix^_i, 

^7 + 1 ^7 ^7 ^7 — 1 

-h^ ^ U2y{x^,yj^i) + hi U2y{x^,yj_^i) 

^-2 ^+2 
+K^;^U2x,y[x.,,yj) - hi^^U2x,y{Xi, yj -l) + kj-'^Ux,2y{Xi,yj) - kj^^Ux,2y{Xi- 
^+2 _ ^_2 ^+2 

— ^ Px {xi,yj)u{xi,yj) — kj 2 Pij "^xixi^yj) 

+kj^^PijUx{xi^i,yj) - hi py{xi,yj)u{xi,yj) 

kf kf 
-hi—^PijUy{xi,yj) + h.,—^pijUy{xi,yj^i) 

-/^.(4,+i - + - 4) + 

(160) +!7i +(7^ -7V^-+A,,+i3,,- 

where 

(161) \A^j\<ckjh^ ^»-i||u||2,oo,a, 

(162) <c/i,fc7'fc^_ih||2,oo,o- 
To simplify the expressions, let 

(163) &.,- = -fc /'+^~^'^ -- 



2 



(164) Cy = -fej ^''^^^ (-/(a;j,?/j+i) +p(xj,?/j+i)uy - w^) 
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^+2 _ ^-2 ^_{_2 ^_2 

—kj — Px{xi,yj)uij — kj - pijdiUij + kj — pijd-^Ui-.ij 

kf - kf kf kf 

(165) -hi py{xi,yj)uij - hi—^pijdfuij + hi—^pijdfu^j^i 

Now, the first correction = (ujj) can be defined as follows 

(166) ul^ = Um+ij = u-o = uIn+i = 0, 
and for £ {1, M} x {1, N}, we have 



(167) (A^Uj ),y + hikjPijUij = / fdxdy + b^j - 6j-ij + cy - c.jj_i + d 



Using inequalities (lEHl), lEU, lESl, (ESI, Ip ^ . lp^ - fTni^ . (fH^ . ltTM| and (fTH2)l combined with the 
estimates obtained in lemma 1301 yields the following 0{h) improvement 

Theorem 3.4. // the solution u of the equation 1^13'^ belonging to C^(f2) and p G C'^{Vl). Let 
be the basic finite volume solution of boundary condition \136\} and discrete equation \13T^ . 
Then the finite volume approximations = {ujj ) defined by the boundary condition i|J6'6]) and the 
discrete equation satisfies the following estimates 

(168) l|e'l|i,r <c//^'||ii||4,oo.a, 



(169) \^k,h,+ .iel^f] <ce 



4.oo,nj 



(170) lle^lU^ <c/i'lkll4.oo.o, 

where ej, = u{xi, yj) — uj, for {i,j) £ {1, M} x {1, N} and vanishes elsewhere. 
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4. Numerical Tests 



4.1. In one Dimensional Space. In this subsection, we give two numerical tests justifying our 
theoretical results in one dimensional case. We mean by uniform mesh that so-called modified finite 
volume scheme and satisfying h^^i = h, h^^^ ~ (see remark 2.5 in T* ) , and by cell-centered 

mesh that satisfying hi = I ^' * ^^aa^ ^^^^ ^ ■ Note that, the first correction computed 

^ "2" ' ^ OClCl 

in case of uniform mesh is the second one, because — /i^ = . 

To show the convergence orders of the first correction and the basic finite volume solution, we 
compute the ratio 

^^^^^ ^ \og{e{h)) - log(e(/io)) 
\og{h) - log(ft.o) 

where ho is the initial value of h in each numerical test and e{K) is the error corresponding to h. 
In the uniform mesh, we use the rule 



ratio = 



log(e(l/2fc+i))-log(e(l/2fc)) 
log2 



4.1.1. First Test. We consider the homogeneous equation (I) : —u^x = f where u{x) = sin{nx) 
and f{x) — t:"^ siniiTx). 

TABLE 1. The convergence orders of the first correction and the basic solution in i^-norm in 

uniform mesh. 



h 


correction 


basic solution 


order 


error //i* 


order 


error //i^ 


1/32 




0.0359 




0.0298e+04 


1/64 


4.0004 


0.0359 


2.0003 


0.1191C+04 


1/128 


4.0000 


0.0359 


2.0001 


0.4764e+04 


1/256 


4.0032 


0.0390 


2.0000 


1.9057e-l-04 



TABLE 2. The convergence orders of the first correction and the basic solution in i^Q-norm in 

uniform mesh. 



h 


correction 


basic solution 


order 


error //i* 


order 


error //i'' 


1/32 




0.1127 




0.0935e+04 


1/64 


3.9999 


0.1127 


1.9999 


0.3742C+04 


1/128 


3.9999 


0.1129 


2.0000 


1.4967C+04 


1/256 


4.0032 


0.1224 


2.0000 


5.9869e-l-04 



TABLE 3. The convergence orders of the first correction and the basic solution in L'^-norm in 

cell-centered mesh. 



h 


correction 


basic solution 


order 


error //i'^ 


order 


error /K'^ 


4/149 




0.2474 




0.3971 


4/599 


1.9903 


0.2507 


1.9986 


0.3978 


4/2999 


1.9943 


0.2516 


1.9991 


0.3981 


4/14999 


1.9961 


0.2473 


1.9994 


0.3968 
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TABLE 4. The convergence orders of the first correction and the basic solution in Hq-hotih in 

cell-centered mesh. 



h 


correction 


basic solution 


order 


error //i"^ 


order 


error /h'' 


4/149 




0.5856 




0.0281e+03 


4/599 


1.9929 


0.5914 


1.0000 


0.1131e+03 


4/2999 


1.9958 


0.5930 


1.0000 


0.5664e+03 


4/14999 


1.9971 


0.5828 


1.0000 


2.8329e+03 



TABLE 5. Comparaison between the accuracy of the first correction that uses first variant and 
the one using second variant in Hq and L^-norms in cell-centered mesh. 



h 


first variant 


second variant 


i^-norm 


i?Q-norm 


i^-norm 


ilg-norm 


4/149 


1.7827e-04 


4.2204e-04 


7.9569e-05 


6.6012e-04 


4/599 


1.1180e-05 


2.6374e-05 


5.4358e-06 


4.1591e-05 


4/2999 


4.4758e-07 


1.0549e-06 


2.2254e-07 


1.6670e-06 


4/14999 


1.7591e-08 


4.1450e-08 


8.6276e-09 


6.6477e-08 



4.1.2. Second Test. In case of cell-centered mest, we saw that the convergence of the first cor- 
rection is the same one of the basic solution in L^-norm for the model (I). We present here an 
example of the mesh where the convergence order of the first correction improves really that one 
of the basic solution in both Hq and L^-norms for the model (I). 
We consider h^_^_l = h and x^_^_l = for all i = 1, N - 1. 

TABLE 6. The convergence orders of the first correction and the basic solution in L^-norm . 



h 


correction 


basic solution 


order 


error /h'' 


order 


error //i"^ 


l/2« 




0.3701 




20.2232 


1/2^ 


1.9668 


0.3788 


1.0000 


40.4472 


l/2l" 


1.9633 


0.3832 


1.0000 


80.8948 


1/2" 


1.9915 


0.3855 


1.0000 


161.7899 


1/2^" 


2.0073 


0.3835 


1.0000 


323.5799 


1/2^^ 


1.9864 


0.3871 


1.0000 


647.1598 



TABLE 7. The convergence orders of the first correction and the basic solution in ifg-norm. 



h 


correction 


basic solution 


order 


error /ft,^ 


order 


error //i^ 


l/2« 




1.4588 




129.5822 


l/2« 


2.1180 


1.3443 


0.9999 


259.1881 


1/2^" 


2.0675 


1.2828 


1.0000 


518.3882 


1/2" 


2.0363 


1.2509 


1.0000 


1.0368e+03 


1/2^" 


2.0301 


1.2251 


1.0000 


2.0736e-F03 


1/2^^ 


1.9985 


1.2263 


1.0000 


4.1471e-F03 
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4.1.3. Third Test. We consider the homogeneous equation (II) : —u^x + + u = / where 
u{x) = sin{TTx) and f{x) = (tt^ + 1) sin(7ra:) + 7rcos(7ra;). 

TABLE 8. The convergence orders of the first correction and the basic solution in L^-norm in 

cell-centered mesh. 



h 


correction 


basic solution 


order 


error /K^ 


order 


error //i"^ 


2/75 




0.1415 




0.0099e+03 


1/150 


2.0074 


0.1401 


0.9796 


0.0407e+03 


1/750 


2.0043 


0.1398 


0.9880 


0.2051e+03 


1/3750 


2.0020 


0.1402 


0.9918 


1.0270e+03 



TABLE 9. The convergence orders of the first correction and the basic solution in i/Q-norm in 

cell-centered mesh. 



h 


correction 


basic solution 


order 


error /K'^ 


order 


error /h''' 


2/75 




0.4441 




0.0377e+03 


1/150 


1.9936 


0.4480 


0.9887 


0.1530e+03 


1/750 


1.9963 


0.4491 


0.9934 


0.7681e+03 


1/3750 


1.9970 


0.4502 


0.9955 


3.8438e-F03 



4.2. In Two Dimensional Space. In this subsection, we present two tests justifying our results 
of two dimensional space. The ratios are computed by using the first formula of ratio in subsection 
5.1 . 



xy{\ — x){\ — y), then u is the solution of H72|l . where 

h, i is even, 



4.2.1. First Test. We consider here u 

f — 2y{l — y) + 2x{l — x). The mesh considered here is such that hi = -l z^'" T "^TV' and 

^ "2" ) ^ IS 00.0. 

kj — 3/1/2, with (xi, yj) is in the center of Kij. 

TABLE 10. Convergence orders of the first correction and basic solution in /fg-norm. 



h 


correction 


basic solution 


order 


error /h'' 


order 


error //?/ 


0.02500 




0.198028 




2.51479 


0.01667 


1.99847 


0.198151 


1.02547 


3.73344 


0.01250 


1.99906 


0.198157 


1.02262 


4.95135 


0.01000 


1.99983 


0.198139 


1.02068 


6.16896 


0.00500 


1.99992 


0.198053 


1.01593 


12.2557 


0.00333 


2.00005 


0.198007 


1.01385 


18.3418 



TABLE 11. Convergence orders of the first correction and the basic solution in L^-norm. 
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ri 


correction 


basic solution 


order 


error /h^ 


order 


error /h'' 


0.02500 




0.0439562 




0.07314 


0.01667 


1.94335 


0.0445737 


1.99851 


0.07318 


0.01250 


1.95651 


0.0448951 


1.99882 


0.07320 


0.01000 


1.96233 


0.045092 


1.99900 


0.07321 


0.00500 


1.97302 


0.0454949 


1.99935 


0.07322 


0.00333 


1.97708 


0.0456318 


1.99947 


0.07322 



4.2.2. Second Test. In this test we choose u{x, y) — sin(7ra;) sin(7r?/), then u is the solution of the 
problem H72|l . where f(x,y) — 27r^ sin(7rx) sin(7r?/). The mesh considered here is such that {xi,yj) 

h, i is even, ^ ^ ( h,j is even, 
I , i is odd ^ \ is odd 

TABLE 12. The convergence orders of the first correction and the basic solution in i/g-norm. 



h 


correction 


basic solution 


order 


error /h'' 


order 


error //i"^ 


0.033333 




2.80465 




34.946 


0.016667 


2.21288 


2.4199 


1.0016 


69.815 


0.008333 


2.17818 


2.1908 


1.0010 


139.59 


0.004166 


2.14753 


2.06364 


1.0007 


279.16 



TABLE 13. The convergence orders of the first correction and the basic solution in i^-norm. 



h 


correction 


basic solution 


order 


error /h'' 


order 


error /h'' 


0.033333 




0.14632 




0.7812 


0.016667 


2.15317 


0.131583 


2.00058 


0.7809 


0.008333 


2.05991 


0.134658 


2.00036 


0.7802 


0.004166 


2.02539 


0.138795 


2.00025 


0.7802 



4.3. Some Comments about the Numerical Results. 

(1) In Table 1 and Table 2, numerical results show that on uniform mesh and for the model 
(I), we can gain an 0{h^)- improvement in both i?Q-norm and L^-norm by the first cor- 
rection. 

(2) In Table 3, numerical results show that for the model (I), we do not have an improvement 
in i^-norm in the first correction when the mesh cell-centered. Furthermore, the cocfficcnts 
of the error in the first correction are better than of those of the basic solution. To improve 
the order in L^-norm, we compute the second correction. 

(3) In Table 4, numerical results show that for the model (I), we gain an 0(/i)-improvement 
by the first correction in iJg-norm. 

(4) In Table 5, numerical results show that the accuracy of the error in the correction defined 
by second variant (see subsection 2.3.1) is batter than that of 117|l in L^-norm and contrary 
in iJg-norm. 

(5) In Table 6 and Table 7, numerical results show that the convergence of the first correction 
improves that of the basic solution in both Hq and norms for (I) . 



is in the center of Kij , hi 
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This implies that, on arbitrarily admissible mesh, the first correction improves the basic 
solution in and norms. 

(6) In Table 8 and Table 9, numerical results show that for the model (II), we gain an 
0(/i)-improvement in both i^-norm and iJg-norm by the first correction for cell-centered 
mesh. 

(7) In Table 10 and Table 12, numerical results show that we gain an 0(/i)-improvement in 
i/p-norm by the first correction. 

(8) In Table 11 and Table 13, numerical results show that, the convergence order of the first 
correction is the same one as of the basic solution in L^-norm, but the errors in the first 
correction are better than of those of the basic solution. 

Remark 4.1. . The idea used here is used by the authors to apply the defect correction technique 
in finite element method with non uniform mesh 
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